Pauli Spin Blockade in the Presence of Strong Spin-Orbit Coupling 
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We study electron transport in a double quantum dot in the Pauli spin blockade regime, in 
the presence of strong spin-orbit coupling. The effect of spin-orbit coupling is incorporated into a 
modified interdot tunnel coupling. We elucidate the role of the external magnetic field, the nuclear 
fields in the dots, and spin relaxation. We find qualitative agreement with experimental observations, 
and we propose a way to extend the range of magnetic fields in which blockade can be observed. 



Blockade phenomena, whereby strong interactions be- 
tween single particles affect the global transport or exci- 
tation properties of a system, are widely used to control 
and detect quantum states of single particles. In single 
electron transistors, the electrostatic interaction between 
electrons can block the current flow [l| , thereby enabling 
precise control over the number of charges on the tran- 
sistor Q. In semiconductor quantum dots, the Pauli ex- 
clusion principle can lead to a spin-selective blockade Q , 
which has proven to be a powerful tool for read-out of the 
spin degree of freedom of single electrons 0, @, 0, 0] ■ 

In this spin blockade regime, a double quantum dot 
is tuned such that current involves the transport cycle 
(0, 1) -> (1, 1) -> (0, 2) -> (0, 1), (n, m) denoting a charge 
state with n(m) excess electrons in the left(right) dot (see 
Fig. [Ha)). Since the only accessible (0, 2) state is a spin 
singlet, the current is blocked as soon as the system enters 
a (1, 1) triplet state (Fig.QJb)): transport is then due to 
spin relaxation processes, possibly including interaction 
with the nuclear fields [j| . This blockade has been used in 
GaAs quantum dots to detect coherent rotations of single 
electron spins 0, 0], coherent rotations of two-electron 
spin states J(|, and mixing of two-electron spin states 
due to hyperfine interaction with nuclear spins 0, Q • 

Motivated by a possibly large increase of efficiency of 
magnetic and electric control over the spin states 13, 11 1 , 
also quantum dots in host materials with a relatively 
large g-factor and strong spin-orbit interaction are being 
investigated. Very recently, Pauli spin blockade has been 
demonstrated in a double quantum dot defined by top 



gates along an InAs nanowire [12j, |13| . However, as com- 
pared to GaAs, spin blockade in InAs nanowire quantum 
dots seems to be destroyed by the strong spin-orbit cou- 
pling: significant spin blockade has been only observed 
at very small external magnetic fields (< 10 mT 12j). 
An important question is whether there exists a way to 
extend this interval of magnetic fields. To answer that 
question, one first has to understand the physical mech- 
anism behind the lifting of the blockade. 

In this work we study Pauli spin blockade in the pres- 
ence of strong spin-orbit mixing. We show that the only 
way spin-orbit coupling interferes with electron trans- 
port through a double dot is by introducing non-spin- 
conserving tunneling elements between the dots. This 



yields coupling of the (1,1) triplet states to the outgoing 
(0, 2) singlet, thereby lifting the spin blockade. How- 
ever, for sufficiently small external magnetic fields this 
does not happen. If the (1,1) states are not split apart 
by a large Zeeman energy, they will rearrange to one 
coupled, decaying state and three blocked states. When 
the external field Bq is increased, it couples the blocked 
states to the decaying state. As soon as this field in- 
duced decay grows larger than the other escape rates 
(i.e. B^T/t 2 > r rc i, where T is the decay rate of the 
(0, 2) singlet, t the strength of the tunnel coupling, and 
r re i the spin relaxation rate [lj]) the blockade is lifted. 
Therefore, the current exhibits a dip at small fields. 

The presence of two random nuclear fields in the dots 
(of typical magnitude K ~ 1 mT) complicates matters 
since it adds another dimension to the parameter space. 
We distinguish two cases: If the nuclear fields are small 
compared to t 2 /V, they just provide an alternative way 
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FIG. 1: Double quantum dot in the Pauli spin blockade 
regime, (a) The double dot is coupled to two leads. Due to a 
voltage bias, electrons can only run from the left to the right 
lead, (b) Energy diagram assuming spin- conserving interdot 
coupling. The only accessible (0, 2) state is a spin singlet: all 
(1, 1) triplet states are not coupled to the (0, 2) state and the 
current is blocked, (c) Energy levels and transition rates as- 
suming non- spin- conserving interdot coupling. We consider 
the 'high'-field limit and neglect the effects of the nuclear 
fields. Then three of the four (1, 1) states can decay, leaving 
only one spin blockaded state |a). Isotropic spin relaxation 
~ P rc i causes transitions between all (1, 1) states. 
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to escape spin blockade, which may compete with spin 
relaxation. There is still a dip at small magnetic fields, 
and the current and width of the dip are determined by 
the maximum of r re i and K 2 T/t 2 . In the second case, 
K 3> t 2 /T, the current may exhibit either a peak or a 
dip, depending on the strength and orientation of the 
spin-orbit mixing. If there is a peak in this regime, the 
cross-over from dip to peak takes place at K ~ t 2 /T. 

Let us now turn to our model. We describe the relative 
detuning of the (1, 1) states and the (0, 2) states by the 
Hamiltonian H e = — A |So 2 ) (<So 2 |, where |<So 2 ) denotes 
the (0, 2) spin singlet state. The energies of the four (1, 1) 
states are further split by the magnetic fields acting on 

the electron spins, H m = Bo{S z l -\-S z r )-\-K l -Sl+K r -Sr, 

where Sur) is the electron spin operator in the left(right) 
dot (for InAs nanostructures g ~ 7 [16]). We chose the 
z-axis along Bq, and included two randomly oriented ef- 
fective nuclear fields K^. R resulting from the hyperfinc 
coupling of the electron spin in each dot to N nuclear 
spins (in InAs quantum dots N ~ 10 5 [l6(, yielding a 
typical magnitude K cx 1/s/N ~ 0.6 peV). We treat 
the nuclear fields classically, disregarding feedback of the 
electron spin dynamics which could lead to dynamical 
nuclear spin polarization [13]. 

Let us now analyze the possible effects of spin-orbit 
coupling, (i) It can mix up the spin and orbital structure 
of the electron states. The resulting states will remain 
Kramers doublets, thus giving no qualitative difference 
with respect to the common spin up and down doublets, 
(ii) The mixing also renormalizes the g-factor that de- 
fines the splitting of the doublets in a magnetic field. 
This, however, is not seen provided we measure Bq in 
units of en ergy , (iii) The coupling also can facilitate spin 
relaxation |l8|, but this is no qualitative change either. 
Some of these aspects have been investigated in (l9| . 

The only place where strong spin-orbit interaction 
leads to a qualitative change is the tunnel coupling be- 
tween the dots. It provides a finite overlap of states dif- 
fering in index of the Kramers doublet (in further dis- 
cussion we refer to this index as to 'spin'), this allow- 
ing for a very compact model incorporating the inter- 
action. The most general non-'spin'-conserving tunnel- 
ing Hamiltonian for two doublet electrons in a left and 

right state reads H t = J2 a ,p { t a0°'L a d 'Rf3 + ^L^}, 
a,j3 being the spin indices, Gum an d o>l(r) the elec- 
tron creation and annihilation operators in the left (right) 
state, and t LR coupling matrices. We impose condi- 
tions of hcrmiticity and time-reversibility on H t and 
concentrate on the matrix elements between the (1,1) 
states and IS02) in our double dot setup. In the con- 
venient basis of orthonormal unpolarized triplet states 
\T X , V ) = iW/2{|T_) T \T+)}/V5, \T Z ) = |T ), and the 
(1,1) singlet \S), this Hamiltonian reads 

H t = it- \f) (5 02 | + t Q \S) (S 02 \ + h.c, (1) 



with \ f) = {\T X ) , \T y ) , IT*)}. The model therefore adds 
a 3-vector of new coupling parameters, t = {t x ,t y ,t z }, 
to the usual spin-conserving to, the vector being a 'real' 
vector with respect to coordinate transformations. If the 
energy scale of spin-orbit interaction is larger or compa- 
rable to the energy distance between the levels in the dot 
(which is believed to be the case in InAs structures), the 
mixing of the doublet components is of the order of 1. 
Therefore we assume that all four coupling parameters 
are generally of the same order of magnitude to,x,i/,z ~ t. 
As the structure of the localized electron wave functions 
is very much dependent on the nanostructure design and 
its inevitable imperfections, the direction of t is hard to 
predict: therefore we consider arbitrary directions. 

We describe the electron dynamics with an evolution 
equation for the density matrix [if]. Next to the Hamil- 
tonian terms, we complement the equation with (i) the 
rates ~ T describing the decay of |So 2 ) an d the refill to 
a (1, 1) state, and (ii) a small electron spin relaxation 
rate r rc i <C L. The full evolution of the electron density 
matrix then can be written as 

^ = -i[H e + H m + H t , p] + Tp + T rol p. (2) 

Experimentally, the temperature exceeds the Zeeman en- 
ergy [lij , allowing us to assume isotropic spin relaxation: 
each (1,1) state will transit to any of the other (1,1) 
states with a rate L rc i/3. Explicitly, we use T rc \ p — 
-r re ip + ^T xe iJ2 a ,d &dPv%, o-l(r) being the Pauli ma- 
trices in the left (right) dot. 

Motivated by experimental work, we assume that the 
decay rate T of |So 2 ) is by far the largest frequency scale 
in @, i.e. P 3> B , K,t, T re \ (in principle T can be compa- 
rable with the detuning A). Under this assumption, we 
separate the time scales and derive the effective evolution 
equation for the density matrix in the (1,1) subspace 

J = -i[H m + H' u p] - G° ut p + G in p + T rcl p. (3) 
The decay and refill terms are now incorporated into 
GkPmn = 2{<WT„ 2 T 2/ + « 2 T 2m }L/(r 2 + 4A 2 ) 

(4) 

GLV„ = feT Il2 T 2m r/(r 2 + 4A 2 ), 

where T a2 = (a\H t \So2)- The coupling between the dots 
gives also rise to an exchange Hamiltonian (H[)ij — 
4A/(r 2 + AA 2 )T. l2 T 2j , with H' t - G out provided r ~ A. 
The diagonal elements of G out give us the decay rates: 
If we consider \T±) and |To), the three triplet states split 
by an external magnetic field, we find T± = Gj_ u _j: ±± = 
2T(t 2 x + t 2 y )/(T 2 + 4A 2 ) and T ° = G$ 00 = 4rt 2 /(r 2 + 
4A 2 ), all of which are - r so - t 2 /T. 

Let us neglect for a moment the nuclear fields, and 
focus on zero detuning, A = 0. This allows us to grasp 
qualitatively the peculiarities of the spin blockade lifting, 
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determined by competition between the Hamiltonian (~ 
Bq) and dissipative terms (~ i 2 /r,r re i) in Eq. [3] 

At sufficiently large fields, the basis states \Tq) and 
\S) are aligned in energy. The spin-orbit modulated 
tunnel coupling then sets the difference between these 
states, which is best seen in a basis that mixes the 
states, | a) = {t | T ) + i t z \S)}/^t 2 +t 2 z and \/3) = 
{it z \Tq) + to \S)}/y/tQ + t 2 . Now | a) is a blocked state, 
0, while decays with an effective rate 
= AT{tl + t 2 z )/(T 2 + 4A 2 ). In Fig. \Rc) we 
give the energy levels of the five states and all transition 
rates in the limit of 'large' external fields. It is clear that 
the system will spend most time in the state \a). The 
current is determined by the spin-relaxation decay rate 
of this state to any unblocked state, 3r rc i/3 = r re i. Let 
us note that if nj, states out of n states are blocked, such 
a decay produces on average n/rib electrons tunneling to 
the outgoing lead before the system is recaptured in a 
blocked state. Therefore, the current is I/e = 4r re i. 

This picture holds until the decay rates of the three 
non-blocked states become comparable with r re i, which 
takes place at B Q ~ v / T i °Trci- To understand this, let us 
start with considering the opposite limit, B <C ^/r so r rc i. 
In this case all four (1,1) states are almost aligned in 
energy, and the instructive basis to work in is the one 
spanned by a single decaying state \m) = {it ■ \T) + 



*o \S)}/y \i\ 2 +*0' an d three orthonormal states |1), |2) 
and |3) that are not coupled to |So2)- At Bq — three 
of the four states are blocked, and spin relaxation to the 
unblocked state proceeds with a rate r re i/3. A relaxation 
process produces on average n/n& = 4/3 electron trans- 
fers, so that the total current is reduced by a factor of 
9 in comparison with the 'high'-field case, I/e — |r rc i. 
This factor of 9 agrees remarkably well with experimental 
observations (see Fig. 2b in Ref. [l2|). 

We now add a finite external field Bo to this picture. 
Since t is generally not parallel to B , the external field 
will split the states |1), |2) and |3) in energy and mix two 
of them with the decaying state |m). This mixing results 
in an effective decay rate ~ i?g/r so , which may compete 
with the spin relaxation rate r re i. At Bo ~ v / r s °r re i, 
we cross over to the 'high'-field regime described above, 
where only one blocked state is left. Therefore, the cur- 
rent exhibits a dip (suppression by a factor 9) around 
zero field with a width estimated as VL sc T ro i (Fig. [5]). 

Let us now include the effects of the nuclear fields Kl,r 
on a qualitative level. If the fields are small compared to 
the scale t 2 /T, their only relevant effect is to mix the 
states described above. This mixing creates a new pos- 
sibility for decay of the blocked states, characterized by 
a rate ~ i^ 2 /F so . This rate may compete with spin 
relaxation ~ r rc i, and could cause the current to scale 
with Tn and the width of the dip with K. In the oppo- 
site limit, K ;§> i 2 /T, the nuclear fields dominate the en- 
ergy scales and separation of the (1, 1) states at Bo ^ K. 



Then, generally all four states are coupled to \So2) on 
equal footing and the spin blockade is lifted. Qualita- 
tively this situation is similar to that without spin-orbit 
interaction (see Eqs 10-12 in 9]). Without spin-orbit in- 
teraction, an increase of magnetic field leads to blocking 
of two triplet states, resulting in a current peak at zero 
field. With spin-orbit interaction, t x>y still couple the 
split-off triplets to the decaying state. Depending on the 
strength and orientation of t, the current in the limit of 
'high' fields can be either smaller or larger than that at 
Bo = 0, so we expect either peak or dip. If it is a peak, 
the transition from peak to dip is expected at K ~ T so , 
that is, at t ~ V KT. Indeed, such a transition has been 
observed upon varying the magnitude of the tunnel cou- 
pling (Fig. 2 in Ref. [L2J). If we assume that K ~ 1.5 mT 
and associate the level broadening observed (~ 100 fieV) 
with r, we estimate t ~ 8 /ieV, which agrees with the 
range of coupling energies mentioned in [121 ] . 

Let us now support the qualitative arguments given 
above with explicit analytical and numerical solutions. 
The current through the double dot is evaluated as 
I/e = P22T, P22 being the steady-state probability to be 
in | S02), as obtained from solving Eq. [2] We give an an- 
alytical solution for A = 0, neglecting the nuclear fields, 
and expressing the answer in terms of the dimensionless 
parameter t/to = ff. Under these assumptions, we find 
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(5) 



with B c = 2V2(1 + \rf\ 2 )(vl + ^r^toVWT and 
/ max = 4er rc i. The current exhibits a Lorentzian-shaped 
dip (see Fig. [21 compare with Fig. 2b in Ref. [12]). The 
width B c and the limits at low and 'high' fields agree 
with the qualitative estimations given above. 

To include the effect of the two nuclear fields, we com- 
pute steady-state solutions of ([2]) and average over many 
configurations of if l,_r @. In Fig. [3] we present the result- 
ing current versus magnetic field and detuning for three 
different regimes. To produce the plots we turned to 
concrete values of the parameters, setting F = 0.1 meV, 
r rc i = 1 MHz, ff = 0.25 x {1,1,1}. We averaged over 
5000 configurations of Kl r, randomly sampled from a 
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FIG. 2: Current as a function of Bo, at A = 0, and neglecting 
the nuclear fields. Around zero field a dip is observed, its 
width depends on the magnitude and orientation of fj. 
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FIG. 3: The current I = ep22T for (a,b) large, (c,d) intermedi- 
ate, and (e,f) small tunnel coupling. The dip observed around 
zero field (a) disappears when io/Y ~ K (c) and evolves into 
a peak for even smaller tunnel coupling (e). 



normal distribution with a r.m.s. of 0.4 /ieV. 

In Fig.[3]Ja) and (b) we assumed large tunnel coupling, 
tl/Y = 6 fieV so that KY/t% = 0.07 is small. In (a) we 
plot the current at A = 0, while in (b) we plot it ver- 
sus detuning for different fixed Bq. We observe in (a) 
a Lorentzian-like dip in the current at Bq — 0. While 
it looks similar to the plots in Fig. [2 the width is de- 
termined by the nuclear fields since K 3> r re i. The 
curve can be accurately fit with the Lorentzian giving 
B c = 7.4 K and I max = 0.62 K 2 T/t . Fig.^b) illustrates 
the unusual broadening of the resonant peak with respect 
to its natural width determined by T. The width in this 
case scales as ~ t^/K and is determined by competition 
of r so and Tat. These plots qualitatively agree with data 

In Fig. Etc) and (d) 
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presented in Fig. 2b in Ref. 
we present the same plots, for smaller tunnel coupling, 
tl/Y = 0.2 fieV = 0.5 K. We included in plot (c) the 
curves for two random nuclear field configurations: It is 
clear that the current strongly depends on Kl.r, which 
agrees with our expectation that in the regime r rc i < Tjv 
the current I oc Tjv oc K 2 . Remarkably, averaging over 
many configurations smooths the sharp features at small 
Bq (c.f. [§]). Plots (d) exhibit no broadening with re- 
spect to T, in correspondence with Fig. 2a of Ref. [lij ]. 
In Fig. [3te) and (f) we again made the same plots for 
yet smaller tunnel coupling, t^/Y — 2 neV <C K. Since 
the nuclear fields now dominate the splitting of the (1, 1) 
states, we see a peak comparable to the one in Fig. 4 of 
Ref. [9( surmounting a finite background current due to 
spin-orbit decay of the split-off triplets. 

We expect our results to hold for any quantum dot 
system with strong spin-orbit interaction. Indeed, recent 
experiments on quantum dots in carbon nanotubes in the 
spin blockade regime [HI display the very same specific 



features, as e.g. a zero-field dip in the current. 

Now that we understand the origin of the lifting of spin 
blockade, we also propose a way to extend the blockade 
region. If one would have a freely rotatable magnet as 
source of the field Bo, one would observe a large increase 
in width of the blockade region as soon as Bq and t are 
parallel. One can understand this as follows. If t effec- 
tively points along the z-direction, t x and t y and thus 
r± are zero: the states \T±) are blocked (see Fig. [2J . As 
\T±) are eigenstates of the field B , this blockade could 
persist up to arbitrarily high fields. Since \Tq) and \S) 
are rotated into \a) and \j3), current will then scale in 
general with the anti-parallel component of spin instead 
of only the spin singlet. 

To conclude, we presented a model to study electron 
transport in the Pauli spin blockade regime in the pres- 
ence of strong spin-orbit interaction. It reproduces all 
features observed in experiment, such as lifting of the 
spin blockade at high external fields or at low interdot 
tunnel coupling. We explain the mechanisms involved 
and identify all relevant energy scales. We also propose 
a simple way to extend the region of spin blockade. 

We acknowledge fruitful discussions with A. Pfund, S. 
Nadj-Perge, S. Frolov, and K. Ensslin. This work is part 
of the research program of the Stichting FOM. 
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